Impact of initial fluctuations and nuclear deformations in isobar collisions* 
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Relativistic isobar (2¢Ru+{$Ru and $871r+98Zr) collisions have revealed intricate differences in their nuclear 
size and shape, inspiring unconventional studies of nuclear structure using relativistic heavy ion collisions. In 
this study, we investigate the relative differences in the mean multiplicity (Riy,,,)) and the second- (R..) and 
third-order eccentricity (Re, ) between isobar collisions using initial state Glauber models. It is found that initial 
fluctuations and nuclear deformations have negligible effects on Riv.) in most central collisions, while both 
are important for the Re, and R.,, the degree of which is sensitive to the underlying nucleonic or sub-nucleonic 
degree of freedom. These features, compared to real data, may probe the particle production mechanism and 


the physics underlying nuclear structure. 
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I. INTRODUCTION 


The recent isobar data [1, 2] stimulated a wide interest from 
the relativistic heavy ion community in the physics of nuclear 
structures of isobars (nuclei with equal number of nucleons 
but different numbers of protons and neutrons) [3-12]. The 
isobar collisions of {$Ru+9$Ru and 96Zr+96Zr were initially 
proposed to search for the chiral magnetic effect (CME) [13- 
16]. The measurements in isobar collisions at \/5,,, = 200 
GeV by the STAR Collaborations showed sizable differences 
in the multiplicity distribution, elliptic flow, and triangular 
flow [1], indicating differences in the density profiles of the 
colliding nuclei [3, 6]. Nuclear density is usually obtained 
from theory calculations, such as energy density functional 
theory (DFT) [17-20], incorporating information from ex- 
perimental measurements. In relativistic heavy ion collision 
simulations, it is common to use the parameterized nuclear 
density in the Woods-Saxon formula, 


PO 
exp [(r — R(1 + B2Y 20 + B3Y2 +...) /a] +1’ 
() 
with the radius parameter R, the diffuseness parameter a, and 
the quadrupole (octupole) deformation parameter (2 (83). A 
large a value in °°Zr is suggested by the STAR data on multi- 
plicity distribution and elliptic flow ratios in non-central col- 
lisions [1], consistent with a halo-type neutron skin thickness 
as predicted by the DFT calculations [21]. The flow ratios 
observed in the most central collisions suggest that the °°Ru 
has a larger 32, while °°Zr has a larger (3 [6]. 
Much progress has been made in quantitatively assessing 
the differences in size and shape of isobars by using relativis- 
tic isobar collisions [5, 7-9, 12, 22-26]. These differences 
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can place stringent constraints on key parameters in theo- 
retical modeling of nuclear interactions, which is the back- 
bone of all theoretical calculations of nuclear structure in- 
cluding DFT. One of such parameters is the slope parameter 
L(p) of the nuclear symmetry energy as a function of nu- 
clear density, which cannot yet be well determined by low- 
energy nuclear experiments [27-35]. In addition, nuclear de- 
formations are largely uncertain from nuclear structure calcu- 
lations [15, 36, 37]. On the other hand, the STAR experiment 
has collected approximately 2 billion minimum-bias events 
for each species of {$Ru+{$Ru and 96Zr+96Zr collisions at 

Syn = 200 GeV, with good control of systematic uncer- 
tainties [1]. These provide an unique opportunity to study the 
density slope parameter of the symmetry energy and the de- 
formation parameters of the isobars. Previous studies suggest 
that ratio observables in isobar collisions are less sensitive to 
final-state interactions [5, 9, 21]. In this work, we investigate 
the effect of initial-state fluctuations on those ratio observ- 
ables in relativistic isobar collisions. 

We use the Optical Glauber model and Monte Carlo 
Glauber model [38—40] to investigate the effect of nuclear de- 
formations and initial fluctuations on the relative differences 


Ri Na) Rea» and Re, in relativistic isobar collisions, where 


Rx = 9 XRuRu _ X ZrZr , 
XRuRu =F X ZrZr 


(2) 


The Glauber model is widely used to determine the central- 
ity in experiments, including the blind analysis of the STAR 
isobar data [1]. Due to symmetry, the elliptic flow v2 with 
impact parameter b = 0 and the triangular flow v3 can not be 
generated in spherical nucleus-nucleus collisions with smooth 
initial conditions in Optical Glauber simulation. The colli- 
sion configuration of deformed nuclei will contribute to spa- 
tial anisotropy, while event-by-event fluctuation contributions 
are also important. We further investigate the effect of initial 
fluctuations by comparing the results under the scenarios of 
nucleon participants and quark participants. 

The rest of the paper is organized as follows. Section II 
gives a brief description of the models and the nuclear density 
profiles used in this work. Section II discusses the effects of 


nuclear deformation and initial fluctuations on Ry Nen)? Fess 
and R.,. A summary is given in Sect. IV. 


Il. GLAUBER MODEL AND ISOBAR NUCLEAR 
DENSITIES 


Particle production in relativistic heavy ion collisions is 
usually parameterized by the two-component model, 


Nen = Npp|(1 — £)Npart /2 + &tNeou]; (3) 


with the number of participants Npart and number of binary 
collisions N.o1 obtained from the Glauber model [39]. Here 
the z is the relative contribution of hard processes and npp = 
Nen in nucleon-nucleon collisions. For the case of x = 0, 
the particle production is independent of hard processes and 
it can alternatively be described by a variant of the Glauber 
model, the so-called Trento model [41, 42], 


p p 1/p 
Nax f s(ra)dira x f (LE Ta) dry 
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(4) 
with p = 1 (we note that previous studies favor p ~ 0 [42]). 
Here s is the entropy density, r} is the transverse radius, and 
T'4(p) is the reduced thickness function for the colliding nu- 
clei A(B) with f (Ta + Tp)d?r1. = Npart. In this work, we 
focus on the centrality dependence of the ratio observables in 
isobar collisions, so the parameter npp and the normalization 
factor hidden in the above equation are not relevant to our 
study. 

The elliptic flow (v2), triangular flow (v3) are the Fourier 
coefficients of the azimuthal angle (@) distribution vn = 
(cos [n (@ — Y,,)]), where (...) denotes the average over an 
event and W,, is the plane angle in momentum phase. The 
large anisotropic flow can be successfully described by rela- 
tivistic hydrodynamics with a large spatial eccentricity [43], 


: n inp 2 
ene?” = frie w(r1)d rL , (5) 
friw(ri dry 


Here w(r 1 ) is the weight factor, often taken to be the entropy 
density s or the energy density e. ®,, is the nth order plane 
angle in the configuration space. The anisotropic flow is ap- 
proximately proportional to the initial eccentricity for small 
amplitude [44]. In this work, we use the participant density 
to calculate €2 and €3, however, our main conclusions do not 
change if alternative w is used. 

For spherical nucleus-nucleus collisions, the multiplicity 
and eccentricities can be calculated analytically in the Optical 
Glauber model for a given impact parameter b. For deformed 
nucleus-nucleus collisions, we calculate the ratio observables 
at a given shape orientation with a specific Euler rotation and 
then take the average over all orientations 


(X(b)) = as J Ordea Ø 


Here X is the observable of interest and Q 4(p) is the orien- 
tation solid angle for colliding nuclei A(B). 


TABLE 1. Woods-Saxon parameterizations of the °°Ru and °°Zr 
nuclear density distributions. The quoted values for Ro and a are in 
fm and that for po is in 1/fm*. 

TRY 67, 
Po R a Bo Po R a B3 
0.159 5.093 0.487 0.00 [0.163 5.022 0.538 0.00 
0.159 5.093 0.471 0.16 [0.163 5.021 0.517 0.20 


The event-by-event fluctuations are not included in the Op- 
tical Glauber model [39], and such an effect is crucial for the 
anisotropic flow observables, especially the odd-order ones. 
We, therefore, use the Monte Carlo Glauber model to include 
the fluctuation effect, where the finite number of nucleons (96 
for the Ru and Zr isobars) are sampled in each colliding nu- 
cleus. The event-by-event fluctuations depend on the number 
of nucleons in the colliding nuclei, so these effects are ex- 
pected to be larger in the isobar collisions than in Au+Au 
collisions. After sampling the nucleons with the appropriate 
nuclear density, the multiplicity is calculated based on Eq. (3), 
and the eccentricity is calculated from the transverse position 
of those finite participating nucleons. Here the participants 
and binary collisions are determined from nucleon-nucleon 
interactions which happens when the distance between the 
two nucleons, each from the incoming and outgoing nuclei, 
is less than \/oy n/r [39], with the inelastic cross section 
set to onn = 42 mb for \/5,, = 200 GeV. Poisson-type 
multiplicity fluctuations are applied in the Monte Carlo sim- 
ulations. 

The density profiles of 9°Ru and °°Zr nuclei are required 
in the Glauber simulations. In this study, the parameters for 
the Woods-Saxon density distributions of °°Ru and °°Zr are 
taken from Ref. [22, 45] (see Tab. 1), which are extracted 
from the isobar densities obtained by DFT calculations with 
the symmetry energy FEsym (pe) = 26.65 MeV and its density 
slope parameter L(p.) = 47.3 MeV at subsaturation cross 
density pe ~ 0.11 fm~? [20]. Nuclear deformation is not 
included in the DFT calculations due to its strong model de- 
pendence. Following our previous work [22, 45], to introduce 
the nuclear deformation, the WS parameters R and a for the 
given (2 (or 33) are calculated to match the volume and root- 
mean-square (RMS) of the corresponding nucleus calculated 
by DFT, keeping the normalization factor po fixed. 

Centrality is determined by the multiplicity distributions in 
the Monte Carlo Glauber model, similar to what was done in 
heavy ion experiments [1]. For the optical Glauber model and 
optical Trento model in spherical nucleus-nucleus collisions, 
the centrality is related to the differential cross-section in im- 
pact parameter [39]. In this study, we use the 16 equal-size 
centrality bins from 0 — 80% centrality range. The results 
obtained with the fixed impact parameter corresponding to 
the centre of each centrality bin are used to mimic the ob- 
servables in such a centrality bin. The impact parameters for 
°6Ru and °©Zr are listed in Tab. 2. The impact parameters are 
larger in 96Zr+96Zr collisions than in {$Ru+{$Ru collisions 
for a given centrality bin because the total cross section for 
6 Zr+I6Zr collisions (o4'2" = 4.551 barn) is larger than that 
for {$Ru+9§Ru collisions (o R"R" = 4.369 barn) based on the 


TABLE 2. The fixed impact parameter used to mimic the 16 equal-size centrality bins from 0 — 80% centrality range. 


Centrality 2.5 TS 12.5 17.5 22.5 275 32.5 
(%) 


37.5 42.5 47.5 52.5 57.5 62.5 


67.5 72.5 77.5 


bru (fm) 
bzr (fm) 


1.865 3.230 4.170 4.933 5.594 6.184 6.723 7.222 7.688 8.128 8.545 8.943 9.323 9.689 10.042 10.384 
1.903 3.296 4.255 5.035 5.709 6.312 6.862 7.370 7.846 8.295 8.721 9.127 9.515 9.888 10.249 10.600 


nuclear density parameters listed in Tab. 1. For simplicity, the 
impact parameter list is also used to determine the centrality 
for deformed nucleus-nucleus collisions in Optical Glauber 
simulations. We note that determining centrality in the de- 
formed case is in principle more complicated, but this does 
not affect our conclusion (see discussion in the next section). 

In order to achieve high precision of calculations in our 
study, the GPU parallel computing technology [46] is used in 
our Glauber model simulations. 


HI. RESULTS AND DISCUSSIONS 


The Finn) obtained from the Optical Glauber model with 
x = Oand x = 0.15, as well as those from the Optical Trento 
model with p = 0 and p = —1, are shown in Fig 1. The 
case of the Optical Glauber simulation with x = 0 is nothing 
but the Optical Trento simulation with p = 1, and the multi- 
plicity, in this case, is only proportional to the number of par- 
ticipants. The Rin) distribution depends on the model pa- 
rameters x and p. The magnitudes increase with x and p over 
the whole centrality range. For the most central collisions, 
the Riya) increases with x in Optical Glauber simulations, 
but is independent of p in Optical Trento simulations. This is 
shown in the top 5% centrality range in Fig. 1(b). The mag- 
nitudes of Ry.) converge to the same value from Optical 
Trento simulations with different values of p, but the trends 
in centrality differ. 

The parameters x and p are usually constrained by the in- 
dividual multiplicity distributions (in this case, of isobar col- 
lisions or their mean). The Rya) in the most central colli- 
sions has been proposed to probe the neutron skin and sym- 
metry energy. However, the particle production mechanisms 
will introduce uncertainties as shown in previous study [5]. 
The results shown in Fig. | imply that the model uncertainty 
of Rina) 1s negligible if the particle production follows the 
Trento model. They also suggest that the precision measure- 
ments of the mean multiplicity distribution of ?7Ru+{$Ru and 
36Zr+36Zr collisions and their difference may offer insights 
on the particle production mechanisms in relativistic heavy 
ion collisions. 

To study effects from event-by-event fluctuations, we 
present the Riy) distributions from Monte Carlo Glauber 
simulations with x = 0.15 and x = 0 in Fig. 2. The apparent 
large fluctuation effect (i.e. difference in Riy,,,) between the 
Optical and the Monte Carlo Glauber) is because of the dif- 
ferent total cross sections predicted by these two models [39]. 
For instance, the impact parameter cut would be smaller for a 
given centrality when a smaller total cross-section is used. 
Based on the nuclear density profiles in our study, the to- 
tal cross sections for {?Ru+{$Ru and 9$Zr+96Zr collisions 


are oy URY — 4.906 barn and oy O% = 4.348 barn, 
which are 4.46% and 3.73% smaller than those in the Optical 
cases, respectively. To remove this trivial difference, we cal- 
culate the optical model using the total cross sections from the 
Monte Carlo Glauber. The results are shown in Fig. 2; they 
become similar now. Our results indicate that both the data of 
centrality dependent R; y.) and the model prediction of total 
cross section are important for the model to data comparison 
in relativistic isobar collisions [1]. 

The ?$Ru and 9$Zr nuclei are deformed [47, 48], as also 
suggested by recent studies [1, 6]. The quadrupole deforma- 
tion parameter for °°Ru is 32,2, = 0.16 and is negligible for 
°67r, The octupole deformation parameter is negligible for 
Ru and is 33,7, = 0.20 for Zr. Figure 3 shows the Rv.) 
from deformed isobar collisions, where the volume and RMS 
radius of the colliding nuclei are constrained to match those 
of spherical nuclei [22, 45]. For both the Optical and Monte 
Carlo calculations, the results from deformed isobar colli- 
sions are similar to those from spherical ones. Here we have 
used the same impact parameters listed in Tab 2 for deformed 
isobar collisions. As mentioned above, the centrality defini- 
tion in the deformed case may be biased by multiplicity vari- 
ations at a given impact parameter. However, the centrality 
is well defined in the Monte Carlo Glauber simulation, and 
the consistency in the corresponding results between spher- 
ical and deformed cases indicates that our centrality treat- 
ment in the Optical Glauber simulation is appropriate. The 
RNa) is considerably larger in the Optical case because of 
the aforementioned larger total cross section. The differences 
are ~ 10% between with and without nuclear deformation. 
Based on previous study [5], this would result in an uncer- 
tainty on the extracted slope parameter of symmetry energy 
of about 5 MeV. Our results indicate that, with appropriate 
magnitudes of deformation parameters, the effect of nuclear 
deformations on the uncertainties of F y) distributions are 
small. 

This is not the case for the anisotropic flow observables, 
where previous studies suggest that both nuclear deformation 
and event-by-event fluctuations are important [39, 49, 50]. 
Apart from nuclear deformation, most of the anisotropic flow 
produced in the most central heavy-ion collisions is due to 
the position fluctuations of nucleons in the colliding nuclei. 
In this paper, we study the effect of nuclear deformations and 
initial fluctuations on geometry eccentricities Re, (n=2,3). 
The results are shown in Fig. 4. Note that, for the Opti- 
cal Glauber simulations, €3 vanishes for the spherical case 
due to symmetry, and for the deformed case the Re, < —1 
(because a > ee | is outside the frame of Fig. 4(b). 
Comparing the results between deformed and spherical nu- 
clei, the general features are similar for both. This is consis- 
tent with the results shown in previous studies [6, 8], where 
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Fig. 1. (Color online). The centrality dependence of R; N.) from Optical Glauber (Opt.Glb) simulations with x = 0.15 and x = 0, and from 
Optical Trento (Opt.Trento) simulations with p = 0 and p = —1 for (a) 0 — 80% centrality range and (b) top 5% centrality. 
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Fig. 2. (Color online). The centrality dependence of R; yN) from 
Optical Glauber (Opt.Glb, solid curves) and Monte Carlo Glauber 
(MC.Glb, symbols) simulations with x = 0.15 and x = 0. The 
dashed curves are the results from Optical Glauber calculations with 


the total cross section from Mont Carlo Glauber simulations. 


the large quadrupole deformation 2 in Ru increases the Re, 
in most central collisions while the large octupole deforma- 
tion 63 suppresses the values in mid-central collisions. 

The magnitudes of Re, show significant differences be- 
tween the Optical and Monte Carlo Glauber models. The Re, 
caused by the nuclear geometry is significantly diluted by the 
event-by-event fluctuations. The reason is rather straightfor- 
ward: event-by-event fluctuations contribute to a large frac- 


tion of the eccentricities in Monte Carlo simulations, and 
these contributions are common between ?¢Ru+{$Ru and 
16Zr+96Zr collisions, suppressing the differences caused by 
nuclear deformation and/or collision geometry. 


In the Monte Carlo Glauber simulation, the participating 
nucleons are usually treated as point-like particles as we did 
in the simulations above. However, nucleons have substruc- 
tures [51-53]. More substructure constituents are expected 
to suppress the fluctuations. We use the quark participant as- 
sumption as it was done in Ref. [51], increasing the density 
by three times and decreasing the inelastic cross-section by a 
factor of 9. The results are also shown in Fig. 4. The Monte 
Carlo Glauber simulations under quark participant scenario 
significantly increase the Re, and decrease R,,, bringing the 
values closer to the Optical limit. We note that despite of the 
cross-section differences, the R,,, calculated in the Optical 
Glauber model are similar to those calculated in the Monte 
Carlo Glauber model with respect to reaction plane instead of 
participant plane [3], which is more related to the elliptic flow 
measurement with respect to zero-degree calorimeter (ZDC) 
instead of time projection chamber (TPC) in experiment [1]. 


There are more sources of initial fluctuations in heavy ion 
collisions, for example, the correlations between the quarks 
in a nucleon, as well as the correlations between the nucleons 
in the colliding nuclei. We find that the minimum distance 
between two nucleons in simulations, as one of the sources 
of nucleon-nucleon correlations, does not change our results. 
We note that the final-state effects, which are not the focus of 
the present study, are also important for the accurate extrac- 
tion of the symmetry energy slope parameter. We postpone 
such investigation to a future work. 
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Fig. 3. (Color online). The centrality dependence of Ryn... with spherical (Sph) and deformed (Def) colliding nuclei, calculated with (a) 
Optical Glauber (Opt.Glb) (b) Monte Carlo (MC.GIb) Glauber models. 
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Fig. 4. (Color online). The centrality dependence of (a) Re, and (b) Rez. The results are calculated with spherical (Sph) and deformed 
(Def) colliding nuclei using Optical Glauber (Opt.Glb) and Monte Carlo Glauber (MC.GIb) models. The results from Monte Carlo Glauber 
model with quark participant assumption (MC.qGIb) are presented as open symbols. From Optical Glauber simulations, €3 vanishes due to 
symmetry in spherical cases, and the Re, < —1 in deformed cases are outside the frame of the plot. 


IV. SUMMARY e RiNa) depends on particle production mechanisms 
and model parameters, while it converges to a single 
value in most central collisions in the Trento model (see 


In this work, the effects of nuclear deformations and initial Fig. 1). 
fluctuations on the ratio observables in relativistic isobar col- 
lisions are studied using the Optical Glauber and Monte Carlo e Initial fluctuations affect the total cross section and thus 
Glauber models. The GPU parallel computing technology is the Rina)» however the effect is negligible for most 
used in our calculations to achieve high precision, especially central collisions (see Fig. 2). 


for the deformed colliding nucleus. Our main findings are as 
follow. e Because of the constraints of the total volumes and the 


measured RMS radii of the isobar nuclei, their defor- 
mations are found to have negligible effects on F Nn) 
(see Fig. 3). 


e The Re, and R,,, driven by the nuclear diffusenesses 
and deformations, are suppressed by the initial finite 
number fluctuations, dependent on the degree of free- 
dom of nucleons or constituent quarks (see Fig. 4). 


Our results indicate unique (in)sensitivities of the ratio ob- 
servables of F yn) and Re, to nuclear deformations and ini- 
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